Hybridized quadrupole-dipole exciton effects in Cu 2 - Organic Heterostructure 
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In the present work we discuss resonant hybridization of the IS quadrupole Wannier-Mott ex- 
citon (WE) in a CuiO quantum well with the Frenkel (FE) dipole exciton in an adjacent layer of 
organic DCM2:CA:PA. The coupling between excitons is due to interaction between the gradient of 
electric field induced by DCM2 Frenkel exciton and the quadrupole moment of the IS transition in 
the cuprous oxide. The specific choice of the organic allows us to use the mechanism of 'solid state 
solvation' (Bulovio^) to dynamically tune the WE and FE into resonance during time w 3.3ns (com- 
parable with the big life time of the WE) of the 'slow' phase of the solvation. The quadrupole-dipole 
hybrid utilizes the big oscillator strength of the FE along with the big lifetime of the quadrupole 
exciton, unlike dipole-dipole hybrid exciton which utilizes big oscillator strength of the FE and big 
radius of the dipole allowed WE. Due to strong spatial dispersion and big mass of the quadrupole 
WE the hybridization is not masked by the kinetic energy or the radiative broadening. The lower 
branch of the hybrid dispersion exhibits a pronounced minimum and may be used in applications. 
Also we investigate and report noticeable change in the coupling due to a induced 'Stark effect' from 
the strong local electric field of the FE. We investigated the fine energy structure of the quantum 
well confined ortho and para excitons in cuprous oxide. 



PACS numbers: 74.78.Fk, 72.80.Le, 78.67.De, 71.35.y, 71.35.Aa 
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INTRODUCTION 



Nanometer-sized organic and inorganic semiconductor 
structures have recently been attracting much attention. 
In these low-dimensional systems, there are pronounced 
quantum confinement effects on the electronic and opti- 
cal properties. Synthesizing composite organic/inorganic 
semiconductors is of major importance not only in the de- 
velopment of novel nano-structure materials for electron- 
ics, optics and transport, but also for basic understanding 
of their size-dependent physical properties. Recently, a 
new type of elementary state which can be generated by 
optical excitation was discussed by AgranovicbA This is 
a hybrid exciton which can be obtained from the resonant 
mixing of Frenkel(FE) and Wannier-Mott (WE) excitons 
in organic-inorganic quantum wells by means of dipole- 
dipole interaction across the interface. Many properties 
of this hybrid were predicted. Other realizations for the 
hybrid have been proposed. Examples are hybrid exci- 
tons in an inorganic semi-conducting quantum dot cov- 
ered by an organic layer— and quantum dot-dendrimer 
system^. The energy of the hybrid exciton as well as the 
Green's function matrix elements for different quantum 
dot-dendrimer systems has been calculated. 

In the model of Agranovich et. al. the decisive role in 
implementing the formation of the hybrid state is played 
by the dipole-dipole coupling between semiconductor WE 
and organic FE. It is assumed that there is no direct wave 
function overlap between the excitons on each side of the 
well. The interaction energy takes the form of P (r)-E (r) 
where P (r) is the polarization field due to the organic 
dipoles of the FE interacting with the electric field E (r) 
in the semiconductor from WE. In the work reported 
here we modified this model to consider the quadrupole 



exciton 'yellow' IS* level in Cu-iO, which is well known 
from many studies on bulk C112O. This immediately puts 
our attention on a different interaction term, which is 
now Qi.jjpEj.k. Here the quadrupole field couples to a 
spatially varying (or k dependent) electric field. This 
coupling is the basis of the present work.— In this paper 
we examine the new hybrid exciton which occurs in a 
Cii20/organic heterostructure. Although the oscillator 
strength of the quadrupole transition is three orders of 
magnitude smaller than the corresponding dipole case we 
show here how this can be compensated by strong spatial 
dispersion of the transition. 

In the following section (JTTJ) , generalizing the work£ on 
dipole-dipole hybrid excitons, formed in adjacent layers 
of organic-inorganic hetero-structures, we will here intro- 
duce a modification. The system we will analyze in this 
paper makes use of dipole forbidden, quadrupole IS exci- 
ton in CuiO (WE) coupled to a suitable organic Frenkel 
exciton (FE). From the fact of strong spatial dispersion 
of the quadrupole transitions we anticipate strong wave 
vector and polarization dependence of the dispersion for 
the hybrid. Concrete results will be given below for real- 
istic configuration and values of parameters for particular 
organic materials. 

To the best of our knowledge, there are no experimen- 
tal studies yet on confinement effects for the cuprous ox- 
ide exciton in a quantum well. Therefore in the third 
([Illjl section we present theoretical examination of con- 
finement effects which are essential due to the rather 
small Bohr radius of 15* exciton; central cell and field 
effects on the exciton dispersion. 

In the fourth pV[) paragraph we discuss the choice of 
some proper organic materials to assure resonance be- 
tween FE and WE in the quantum wells. We are going 
to show that layers of PS:CA:DCM3i give an excellent 
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match to the quadrupole 'yellow' exciton of C1L2O. We 
have chosen this organic compound due to three main 
factors: 

1 . the extremely big oscillator strength of the organic 
Frcnkel exciton formed on DCM2 molecules; 

2. unlike the conventional organic with short emission 
lifetime, the dynamic 'solid state solvation' mech- 
anism discovered in such a compound would allow 
the hybrid to live through the phase of 'slow' sol- 
vation (pa 3 ns) which is comparable to the lifetime 
of the quadrupole exciton; 

3. the ability to tune into the resonance the energy of 
the singlet FE simply by means of changing con- 
centration of the CA. 

The strong local electric field induced by the organic 
Frenkel excitons penetrates into the cuprous oxide layer 
and results in an induced 'confined Stark effect'. Rel- 
ative shift of the electron and hole give rise to induced 
polarization and reduced hybridization effect. This will 
be discussed in detail in the fifth (|V| section. 

Hybridization requires coupling between the two exci- 
tons, and we will estimate the coupling coefficient in the 
sixth (|VI[) section. The coupling between FE and WE in 
the case of quadrupole active 15" WE is due to gradient 
of the field induced by the FE in the DCM2 organic. In 
the dipolc-dipole hybrid exciton one utilizes big oscillator 
strength of the FE and big Bohr radius of the WE, then 
in case of quadrupole-dipole exciton one utilizes big os- 
cillator strength of the FE and long life time of the WE 
along with enhanced spatial dispersion of the coupling 
parameter. 

In the final section (|VII|) we will discuss the specifics 
of the quadrupole-dipole hybridization dispersion and 
briefly propose possible applications of such a hybridiza- 
tion. We look forward to experimental tests of our re- 
sults, by means of new types of high precision spec- 
troscopy which were invented recently^. 



II. THE QUADRUPOLE-DIPOLE HYBRID 
EXCITON CONFIGURATION 

Our proposed configuration of an experiment for ob- 
taining the hybrid exciton is shown in the following dia- 
gram (see FigfTJ). In this simplified model a mono-layer 
of width L w pa size of a unit cell 4.6 A quantum well of 
Cu-iO (gap energy E g — 2.17 eV—) is placed upon a thin 
film of the PS:CA:DCM2 organic (with the gap energy 
much bigger than that of the cuprous oxide). Obliquely 
incident (to assure x component of the exciton wave vec- 
tor) photons of close energies excite Wannier-Mott exci- 
ton in the cuprous oxide in resonance with the Frenkel 
exciton in DCM2. We consider DCM2 exciton as a 2D 
lattice of dipolcs d x = 12D at discrete sites n, placed at 
z pa L w /2. 




Cu 2 PS :DCM2: CA 



FIG. 1: Schematic representation of a possible experimen- 
tal set-up to produce quadrupole-dipole excitons. Here the 
inorganic C1L2O quantum well provides the Wannier-Mott 
15 quadrupole exciton with the binding energy 143 meV 
(for details see section IIII|) . Two pumping photons Tiuj ~ 
Edcmz, Eis generate the hybrid signal from the upper E u 
and lower Ei branches. The DCM2 part of the organic 'solid 
state solute' provides dipole allowed Frenkel exciton; the PS 
host prevents wave function overlapping between organic and 
inorganic excitons; CA under proper concentration allows 
tuning of the excitons into the resonance. Due to comparable 
lifetimes of both types of exciton the system utilizes strong 
spacial dispersion of the quadrupole exciton and big oscillator 
strength of the organic. 

We treat the interaction of cuprous oxide quadrupole 
excitations only with DCM2 organic molecules in the 
quadruple approximation, since at the pumping laser 
frequency only dipole forbidden transitions are allowed. 
Due to the small concentration of the DCM2 molecules 
there is a 'buffer' of PS between Cu 2 and DCM2 so 
that one may neglect the exciton wave function overlap 
and assume perfect 2D invariance of the system in direc- 
tion transverse to growth. Further on we neglect kinetic 
energy of the Wannier exciton. Indeed, as will be shown 
in detail in the next paragraph, due to fluctuations of 
the inorganic quantum well (IQW) width, strong con- 
finement of the exciton occurs^ compared to the kinetic 
energy for small wave vectors which results in hopping 
motion of the exciton between sites of localization. We 
seek the new quadrupole-dipole hybrid eigen-states for 
the upper (u) and lower (I) branches in the usual linear 
combination form: 

|«,k) = A u \F,k}+B u \W,k) , \l,k) = A l \F,k}+Bi\W,k) 

where the weighting coefficients for small k are given by^: 

K,i=Bl,i = 1/2- 

After the usual diagonalization of the coupled WE /FE 
hamiltonian, H = Hwe + Hfe + Hi„ t , the energies of 
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the resulting upper and lower branches are given by: 

where we have introduced the quadrupole-dipole inter- 
action Hamiltonian and corresponding off-diagonal hy- 
bridization parameter: r(k)= (W, k| Hint \F, k) . 

In the subsequent paragraphs we will derive expres- 
sions for the energies of the exciton and necessary con- 
ditions for the resonance hybridization. Also as a main 
result of the article we will derive an analytical expression 
for the hybridization parameter and dispersion. This will 
be followed by quantitative and qualitative comparison of 
systems with dipolc-dipolc hybridization^. 



III. QUASI 2-D EXCITON IN Cu 2 0. 

As far as we are aware there is no experimental data or 
theoretical description for the film cuprous oxide systems; 
contrary to the extensive literature for the bulk case. In 
case of the small size C112O quadrupole exciton one can 
not consider strong confinement effect, which requires the 
width of the IQW much less than the Bohr radius of the 
exciton a^. Even if modern epitaxy methods allows one 
to get a molecular mono-layer for the quantum well thick- 
ness, one would have L w /ab ~ 1/2 which is not enough 
for pure 2D consideration. Two and more mono layers in 
one quantum well will give the case of weak confinement 
and result in much weaker coupling (see below). Thus 
we start from the well described case of bulk "yellow" 
excitons and then estimate the main properties in case 
of strong confinement theoretically. 

CuiO condenses in a cubic structure, where the cop- 
per ions form a face-centered sub lattice, while the oxygen 
ions form a body-centered sub lattice. The arrangement 
of both sub lattices is such that a copper ion is found 
centered between two neighboring oxygen ions (Oh sym- 
metry) with the lattice constant a ss 4.26 A. From the 
lattice structure we now turn to the band structure of this 
crystal^. There is a direct band gap, where the valence 
band is formed by the Cu 3d orbitals and the conduction 
band arises from Cu 4s and (possible) oxygen orbitals. 
When considering the cubic crystal field, the five 3d 
states of the valence band split further into three states of 
type 3 r^ and a twofold 2 r^" level. Taking also spin-orbit 
interaction into account, the state splits further into a 
twofold level 2 T^ and a fourfold degenerate 4 r^~ level. 
The exciton representation is obtained from the direct 
product of electron and hole T ex = T enve i op (&T e (g) T/j. 
For S-excitons : <g) 2 rf ®r£ = 3 T++ 1 T+ The three- 
fold degenerate 3 r^ state and single 1 rJ state are termed 
ortho-exciton and para-exciton respectively. The Para- 
exciton is optically forbidden in bulk. The Ortho-exciton 



^ 3 T 4 ) forbidden from the ground state of sym- 



is dipole 

metry ( 1 r+), because ^Tf \ 3 T^ \ 3 T+) = and couples 
to the light in lowest order via quadruple interaction of 

Cr+| 3 r+ pr+)^o. 



Unlike the dipole operator, the quadrupole operator 
depends on the direction of the light wave vector k rela- 
tive to the lattice and the radiation polarization vector e. 
Because of the k dependence the transition is anisotropic 
even in a cubic crystal. The amplitudes of the ortho- 
exciton quadruple transitions are given by the symmet- 
ric vector product of k and e : ~ e^k, + eykj, i =/= j, see 
for example^. The three components correspond to the 
Cartesian representations: 3 r^. z , 3 T^, 3 . 

The measured^ oscillator strength of the quadrupole 
transition in bulk Cu^O is low / ~ 3.9 • 10~ 9 , which is 
about four orders of magnitude smaller than the value 
found for the dipole transitions of the P-excitons of the 
yellow series. Even though the coupling to the light is 
extremely weak, it can not be disregarded. The binding 
energy is about 153 meV and Bohr radius of the exciton 

is given by a b = 2|la « 7 A. 

Because of the unidirectional confinement in the IQW, 
the exciton is discretized in this direction (z-dircction). 
And so the symmetry group Oh is reduced to D^h- As 
a result the three fold degenerate ortho-exciton 3 r^ is 
split into two fold degenerate 2 r^" and non degenerate 
X T\ ortho-exciton levels, which we are going to refer to as 
"heavy" and "light hole exciton" in analogy with the well 
known case of dipole allowed exciton. Another remark- 
able result of the confinement is that the para-exciton 
l Y\ changes its symmetry to And due to the fact 

that the quadrupole operator 3 r^" also reduces its sym- 
metry to 2 r^" and ^Tjj" the para-exciton is no longer for- 
bidden in the IQW. Note that roughness of the interface 
due to the small width of the IQW leads to further reduc- 
tion of the symmetry to D^h ■ The details of the selection 
rules and dependence upon polarization can be found in 
the Table I, II. 

Now let us consider the effect of confinement on the 
energy spectrum of the quadrupole active exciton. We 
study only the ortho-exciton (the generalization to the 
case of para-exciton is trivial) and neglect spin-dependent 
exchange interaction and related spin-orbit effects. The 
energy of the 15* quadrupole exciton confined in a IQW 
can be written as: 



Hwe = H e (z e )+Hh (zh)- 



p 2 P 2 



2 ^ 2M eJ p 2 + (z e - z h f 



symmetry ( 3 T^ 



where p is the relative electron e - hole h position, p is the 
relative momentum, P is center of mass quasi momentum 
which is to be considered in x, y— plane, (i and M are 
reduced mass and total mass of the exciton respectively. 

The quadrupole exciton confined in the infinite IQW 
has a smaller Bohr radius then in case of bulk and as a 
result one must take into account so called "central cell" 
corrections (CCC) 9 in determining the dispersion (See 
Fig®. 

Aside from the non-parabolicity of the bands 1 ^, the 
electron-hole interaction in this case is the bare Coulomb 
interaction modified by the k dependent dielectric func- 
tion. This is mainly due to virtual LO-phonon assisted 
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FIG. 2: The solid line shows the modified energy of the IS 
exciton as a function of the Bohr radius (A) and variational 
parameter A due to the confinement effect; the dotted line 
shows the same function with the central cell correction not 
taken into account. 



valence electron transition into the highest conduction 
bands ( split by the confinement). Here we have 
to assume that the energy of the LO-phonon modes ex- 
ceeds 87 meV— , when the dielectric constant drops from 
co = 7.5 =F 0.2 to Eoo = 6.46. In this case it can be 
shown^ that for small values of the exciton wave number 
k, e(k) weoo -0.18 (fca) 2 - 

In our model due to the comparable sizes of the exciton 
radius and IQW width one can consider the confinement 
as a small perturbation to the pure 2D exciton energy. To 
explicitly show the perturbation part it is convenient to 
introduce a small variational parameter A ranging from 
1 for pure 2D to 1/2 for pure bulk cases, so that the 
confined exciton energy can be re-written as: 



WE 



H\VE 



H 



WE 



(1) 



Where we separate the analytically solvable part: 

H° WE = H e (z e ) + H h (z h ) + f + ^ ? -— (2) 

and small perturbation parts are due to the weak con- 
finement: 



mass. The fourth term in ((2|) and second term in (|4|) 
yield the free exciton dispersion relation for the center of 
mass motion^: 



21.2 



£k = 



nrk 

2M 



1 - 



2X 2 Ma 2 h 



' k D a b ' 
, 2A / 



M'a 2 



' k D a b ^ 
, 2A / 



(5) 



Where ko is the Debye wave vector, i.e. the radius of 
a sphere with volume equal to that of the first Brillouin 

Kouvalkis et. al£ estimated 



zone and /„ (x) 



i (1+2/ 2 ) 1 " 



the order of magnitude of the constants C e , C'h using k • 
p perturbation theory. Mixing of bands with different 
parities modifies the bare electron and hole masses. In 
addition including the coupling to LO modes gives C' e f=a 
C A wl. 

In our work we are going to use a simplified approach 
instead the k • p perturbation theory and estimate this 
correction from the fact that the exciton is becoming lo- 
calized as the Bohr radius approaches the lattice con- 
stant: lim Ek = 0. This yields the following expression 

a b — >a 

C = 2I2 (tt) /I4 (tt). The correction to the potential en- 
ergy due to the first term in ((4]) is given by its expectation 
value in momentum space. 




Here the re lativ e electron e - hole h motion wave function 
-7 is strongly peaked in momentum space. 



l + {qa b ) 



The last term in Hyy E describes the effective interaction 
between an electron and hole at momentum transfer q 
with the same approximation for the k dependent dielec- 
tric function as above: 



V(q) 



Aire 2 



(9) 



Aire 2 
q 2 £oc 



Aire 2 



0.18. 



11 WE 



£\J P 2 + (z e ~ Z h )' 



If one truncates the IS trial wave function outside the 

(3) 

first Brillouin zone we find the correction due to the small 
Bohr radius of the exciton in the form: 



and 



if 2 

11 WE 



4 2 

pa 
2Ah 2 n' 



p 2 P 2 a 2 
4h 2 M' 



(4) 



. These terms (|3|4|) come from second order corrections in 
the tight-binding model and are due to non-parabolicity 
of the bands with 1/fj,' = C e /m e + Ch/rrih and M' = 
m e /C e + mh/Ch- The second term in ^ couples the 
relative motion of the electron and hole with the motion 
of their center of mass and modifies the total exciton 



0.36A 3 e 2 a 2 [4(^)] 2 

In the above expression the contribution of LO phonons 
to the dielectric function has been employed, where the 
main contribution comes from virtual transitions between 
the higher Tg conduction band and the highest va- 
lence band. Now treating H^e as perturbation in the 
lowest order (neglecting momentum dependence of the 
screening and approximating H e (z e ) + (zh) + 557 by 
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^4- ) one can find the total Wannier-Mott exciton en- 
A b . 

ergy in the IQW as the minimum of the total exciton 
energy with respect to the Bohr radius : 



E 



total 
WE 



(A, Ob) = 



h 2 \ 2 



h 2 a 2 C h 



. 2A ; 



0.36A 3 e 2 a 2 [l' 2 {^t)\ 



(6) 



24m / 2 



' kjyO± 

v 2A 



' koOb 
> 2A 



To find the variational parameter A we use an additional 
restriction determined by requiring that the first order 
energy shift vanishes^ 2 -: 



(2A)^exp(-^) _. kr 
*a = Xe (z) X h (z) ^ b e 



(* A |F^ B |* A )=0 (7) 
(8) 



Where is the unperturbed eigenfunction of H^ VE . 
The envelope functions for the confined electron and hole 
are denoted as Xe and Xh correspondingly and for k = 
given by equation (|14[) . The numerical calculations (see 
fig[2]and figE]) show that for a mono-layer of the cuprous 

o 

oxide A ~ 0.881, ab m a~ 4.6 A, i.e. in this case one can 
consider the IS quadrupole exciton to be Frenkel like 
localized exciton with k dependent oscillator strength. 
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FIG. 3: (Color online) For numerical estimation of the con- 
fined quadrupole Bohr radius as an approximation we took ab 
to be the bulk Bohr radius 5.1 A. So the 'correction' equation 
(0 becomes function of only one parameter A and so does the 
binding energy (J6j . 

Without central cell corrections the value of the bind- 
ing energy is 124 meV (minimum of the dashed curve 
on FigJ3]). With CCC it lowers to 154 meV and it cor- 
responds to the minimum on the lower solid curve of 
Figj2] Now in the weak confinement of the IQW there 
is a restriction on the value of the parameter A given by 
equation ([7]). So the standard binding energy (second 
term in the total energy ([6])) slightly grows with the in- 
creasing parameter A^ 3 -. But the CCC with increasing 
confinement A decreases. This results in binding energy 
for weak confinement bigger than bulk case without CCC 



but slightly smaller then for the bulk case with CCC. 
This unusual behavior is entirely attributed to the con- 
finement dependent CCC. 

The confinement in the IQW increases the overlap of 
the electron and hole wave functions, which in turn in- 
creases the oscillator strength of the 'yellow' transition. 
We are only able to estimate the oscillator strength in 
case of quantum confinement. Wc consider a non ra- 
diative interface exciton, which refers to the states out- 
side the photon cone, k > The exciton propagat- 
ing in the plane are trapped and accompanied by the 
light field which is evanescent in the direction perpen- 
dicular to the interface (if one considers formation of 
polariton modes), i.e., are invisible at macroscopic dis- 
tances from the IQW. In the strictly two-dimensional 
limit, the oscillator strength of the lowest state scales 
f 2D 

Sa b 

scales as 



(A fc ) 



■ . The oscillator strength per unit volume 
4-. The maximum enhancement of the 



oscillator strength in IQW with respect to the bulk case 
is eight times and results in giving f 2D w 8 ■ 3.7 x 10 -9 . 
In general, due to the interaction with the Frenkel exci- 
ton f 2D will have weak dependence on the wave vector. 
Also the exciton resonance broadens due to imperfections 
of the IQW. In our case, the IQW is rather thin, which 
gives rise to interface roughness and thickness fluctua- 
tions. Therefore the exciton mode should show a line 
width of about 1 \ieV (the lifetime ~ 1.7 ns). 



IV. PS:CA:DCM2 AS AN ORGANIC PART OF 
THE HYBRID 

As it was already discussed, for the best manifestation 
of the hybridization effect one has to be able to tune the 
energy of the Frenkel and Wannier-Mott excitons into 
resonance. Also the Frenkel exciton lifetime should be 
comparable to the lifetime of the IS quadrupole exciton, 
otherwise the hybrid would not live long enough to uti- 
lize peculiar properties of the quadrupole part. So we 
decided to use in our model the Frenkel exciton formed 
in a recently reported amorphous organic thin film doped 
with the red laser dye: [2-mcthyl-6-2-(2,3,6,7-tetrahydro- 
1H, 5H - benzo[i,j] - quinolizin - 9 - yl) - cthenyl] - 4H - 
pyran - 4 - ylidene] propane dinitritle (DCM2)i2, which 
has an electric dipole moment ~ WD in the ground 
state. To achieve a red spectral shift of DCM2 into the 
resonance with the IS quadrupole exciton wc propose to 
use a low DCM2 dopant concentration in a two compo- 
nent host consisting of polystycrine (PS) and the polar 
small molecule camorphic anhydride (CA). 

Because DCM2 and CA (dipole moment « 6D) are 
highly polar molecules, and PS is a nearly non-polar (less 
then ID), one can adjust the spectral shift by means of 
adjusting the relative concentration of the DCM2 and 
CA molecules. Increasing the DCM2 concentration one 
increases the strength of the local electric fields present 
in the film. In our case we will keep DCM2 concentra- 
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tion constant and low (0.05%), to avoid overlaping of the 
organic and inorganic excitons, and limiting DCM2 ag- 
gregation effects. At the same time , the dielectric prop- 
erties of the film are modified by changing the concentra- 
ton of CA, which has a large ground state dipole moment 
relative to its molecular weight and is opticaly inactive 
in the relevant region of the DCM2 photoluminiscence. 
The dielectric permitivity increases with increasing CA 
concentration linearly: 

e = 2.44 + 0.13 (CA%) (9) 

In contrast, the index of refraction n m 1.55 of the film is 
nearly constant. The PS provides a transparent , non po- 
lar host matrix. Such a mechanism for the spectral shift 
was termed 'solid state solvation', as the solvation mech- 
anism underlying 'solvatochromism' of organic molecules 
in liquids^. 

The theory of 'solvatochromism' relates the experimen- 
tally observed changes in emission and absorption spectra 
of a solute (DCM2) to the dielectric permitivity of a sol- 
vent (PS:CA). An electronic transition on the solute (due 
to photon absorption) produces a corresponding change 
in the solute charge distribution, which causes the sur- 
rounding solvent molecules to respond to this new field 
apart from the Frank-Condon (FC) shift due to solute 
nuclear reorganization (See Figf5]for details) in two ways: 

1. through electronic cloud reorganization (polariz- 
ability) which is referred to as 'fast' solvation which 
occurs during 0.16 fs judging from the relaxation 
spectrum of the DCM2 (line width w 0.25 eV); 

2. through gross spatial movement due to physical 
translation and rotation ('slow' solvation). Dur- 
ing this phase the radiative recombination from the 
FE is prohibited and allows the FE exciton to live 
3.3 ns. In our proposed sheme , the actual hy- 
bridization with IS quadrupole exciton occurs dur- 
ing this time interval as it is comparable to the life 
time of the quantum confined quadrupole exciton, 
while the lifetime of the DCM2 in vacuum is deter- 
mined by FC effect which is much faster; 

Once the Frenkel exciton energy falls into the non zero 
coupling parameter vicinity with the quadrupole exciton 
(resonance condition) the hybridization occurs. The de- 
tailed dynamics of the hybrid is the subject of our future 
work. Presently we assume that the effective lifetime of 
the hybrid is ~p^- and requires two photons to popu- 
late both branches of the hybrid exciton. We performed 
calculations for the energy of the Frenkel exciton under 
the 'continuum' approximation, in which the surround- 
ing molecules of the solvent are replaced by a continues 
dielectric. The molecules of the solute are described by a 
spherical cavity of radius ap of the DCM2 molecules and 
corresponding charge distribution is reduced to the dom- 
inant dipole moment^. Then the total emission energy 
including solvation effect can be written as: 

E = E + AE solv (10) 



E F , 
\e) 




FIG. 4: Schematic of dynamic tuning of Frenkel and IS 
quadrupole Wannier excitons by means of 'solid state sol- 
vation' (SSS) red shift effect. Following photon absorption 
(1) in DCM2 there are distinct dynamic phases: (2) fast ad- 
justment of the electronic configuration in DCM2 due to in- 
teraction with polar CA molecules within the time frame of 
pico seconds t/; (3) slow self-adjustment of DCM2 and CA 
molecules t s m 3.3 ns; (4) second 'yellow' photon delayed by 
St ~ t s — ns is absorbed by Cu-iO confined 15 quadrupole 
exciton with lifetime t\s ~ 1.7 ns; hybridization occurs when 
detuning between this two types of exciton becomes smaller 
then the k dependent coupling parameter, lifetime of the hy- 
brid is defined by r s and ns; (5) phosphorescence due to the 
hybrid exciton recombination. (6) shows possible recombina- 
tion of non-interacting excitons. 

where Eq is the emission energy in vacuum including the 
FC shift; and E so i v = —-^(jig — ^ e )(A^ e + A op fi g ), here 
we introduced ji e and fi g to be , respectively, solute ex- 
cited and ground states dipolcs, and 

A=*<£zi> | A op = % - 1) 
2^+1 p 2n 2 + l 

associated with 'slow' and 'fast' relaxation processes (de- 
scribed above) respectively. 

The equations (|9|10|) along with experimental fitting 
for —-^-{i_tg — /i e )/i e = 0.57 yields the concentration of 

the CA f=a 22% to correspond to the resonance with the 
confined 15* quadrupole energy 2.05 eV. 

There is one sensitive point we make in our work, 
namely necessity of the second delayed photon. One pho- 
ton does allow hybridization once the energy of the FE 
exciton is close enough. Although only one branch of the 
hybrid is going to be populated. Because we approach 
the resonant energy from above, the upper branch is the 
one to be populated. The lower branch may be popu- 
lated by multi acoustical phonon process. But first we 
do not know how long it takes to populate from the max- 
imum of the upper branch to the minimum of the lower 
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branch. For the temperature of 1.7 K it takes approxi- 
mately 26 of such phonon processes without Stark effect 
(see section E} or 7 of them with the Stark effect. Even 
though such transitions are allowed it may take more 
time than life time of the hybrid. For example the life- 
time of the IS ortho-exciton itself is determined by the 
acoustical phonon assisted transformation to the para- 
exciton which involves change in energy 12 meV. Which 
is comparable to 2T k . 

So the second photon is designed to assure both 
branches are populated regardless of the initial condi- 
tions. When both photons are in resonance with the WE 
exciton hu w E\s, the second photon provides conserva- 
tion of the energy: 

hus + hw = (E 1S + T k ) + (E 1S - r fe ) 



V. INDUCED STARK EFFECT 

In this paragraph we investigate a new effect which 
plays a rather significant role in forming any hybrid ex- 
citon. Let us refer to it as an induced quantum confined 
Stark effect. The main idea is the following. There is 
an evanescent potential in z - direction and correspond- 
ing k - dependent electric field due to the Frenkel or- 
ganic (DCM2) exciton (FE), penetrating into the inor- 
ganic quantum well (see FigQ}. This field induces some 
polarization of the bound electron-hole pairs forming the 
Wannier-Mott inorganic exciton (WE). This polarization 
leads to an effective screening of the electric field and re- 
sults in reducing the coupling. 

We consider the FE as polarization wave^ P (r) con- 
fined to a perfectly 2D organic quantum well (OQW) 
placed at position z' ~ ^f-. Neglecting the constant 

phase of the wave e~' i— ^ the polarization per unit area 
S is given as: 

-tknn / t 

p W-^V(-T 

Now let us we can move from discrete set of FE to a 
continuous distribution and take into account that 

1 

a 2 



a F 



li 



dru 



In this continuous approximation one gets the final result 
for the polarization created by the FE exciton propaga- 
tion along the interface: 



P(r) 



-ikrn 



c.c 



aF\fS V 2 
The potential due to this polarization wave is given by^: 

ip k (r) = G (r, n|| , k) VP = -V,G (r, n„ , k) P, 
The standard Green's function for z < z' is given by: 

4tt e k ( z ~ z ') 



G (z, z , k) 



It gives rise to the potential of the form: 

<^k (r) - J2 ^ + * fc l G ( k ' z '' z ) P * W 

47T 



= (r) = d- 



(11) 



aF 



Here we utilized the fact that the complex conjugate part 
cancels out the imaginary part of the above expression. 
The electric field creates an additional polarization in 
the IQW, which we are going to estimate by its effect 
upon the coupling parameter. Effectively this polariza- 
tion manifests itself in an envelope function equation in 
Ithe QW for the electron and hole. Following a standard 
procedure ( Bastard^) , the equations governing the elec- 
tron (hole) envelope function in the IQW: 



2m e dz 2 

h 2 d 2 

2nih dz 2 



-eF k z^j xi (z) = (E v -E)xt(z) 



Where we used the long wave approximation for the po- 
tential namely taking the effective electric field in the 
form: F k ~ ktfk (z = z'). The total energy shift due to 
the induced field is negligible for small difference in the 
electron and hole masses and small width of the IQW. 
The envelope functions are subject to zero boundary con- 
ditions at both sides of the IQW. Here one comes across 
some difficulties. Although this system has an exact so- 
lution the standard approach with Airy functions Ai and 
Bi will lead us to a rather complicated but exact result: 



Xfc (z) = C lt n tk Ai 



eF k z - Cn.fc 

eF k z - Cn.fc 
Co.fc 



(12) 



In the last expression we introduce the following notation 
for the electron (hole) energy change due to the confine- 
ment and FE induced polarization: 



1/3 



( (eF k h) 
Cn,k = E n -E - eip k (0) 



With E = for the electron and E = —E g for the 
hole (we also omitted indexes e, h for simplicity). The 
normalization constants are connected as: 



Ai 



Co.fs 



(k)+e k 



The energy levels can be found using boundary condi- 
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tions as discrete solutions for: 



-Bi ( eF ^~^A M - eFkL f- 



(13) 



Co.fc 



= 



The constant is defined from normalization condition. 
For the hole one has to change e — > — e. Although (fT3"|) 
and ([12")) exhibit the exact solution, they are rather com- 
plicated for further analytical description of the hybrid 
exciton states. 

There are some approximate ways to treat the problem 
of relative shift of the electron and hole wave functions 
in the induced field. First let us treat the electric field 
Fk using perturbation theory. In this approach we can 
explicitly see the term due to induced polarization. This 
perturbation approach is applicable if the energy differ- 
ence between unperturbed ground and first exited state 
is much bigger then the perturbing potential at the av- 
erage position of the particle z = 0. For the case of a 
mono-layer of C112O it is applicable only for very small 
wave vector region but may be used as a correction in 
the dipolc-dipole hybridization. In the first order there 
will be no change in total energy ( cx F%) but the wave 
functions will be changed to: 



xl (z) 



32F k Ll,em fl . /2ttz 
- I J — \ 7 — sm I "7 — 

L w J 27ir 2 Ti V L w \L W 

(14) 

Because (|14p is still rather complicated, we will consider 
only the lowest energy levels transitions so it is possible 
to consider instead of the infinite IQW, an 'equivalent' 
parabolic profile defined through its lowest energy level 
as 



1 



(jrKf 
2mL„ 



So we are able to consider k, L w values in the region 
where the perturbation theory is not applicable. This 
problem has the exact solution with the same energy as 
was given by the perturbation theory for the lowest tran- 
sition. The envelope function normalized inside the given 
IQW has a Gaussian form: 



(27T) 



1/4 



y/erf(n/V5)L t 



: exp 



(z - z e . h y 

Rl 



(15) 



In (|T5|) the dimension i?o of the parabolic IQW for the 
electrons and holes is taken to be the same: 



Ro — 



m e uj 



7T 



The shifted average electron and hole positions are given 
by: 



Ze,h 



± eFi 



2mujQ 



2w 4 h 2 



In the next paragraph we are going to use these results 
for the envelope function to calculate the coupling pa- 
rameter. 



VI. THE COUPLING PARAMETER AND 
DISPERSION 

The energy of interaction between the organic dipole 
and inorganic 15* quadrupolar excitons can be written 

a all. 



Ha, 



1 



2^ Qof* 



a (3 



d 

dx a 



Fa 



= -g&*,/3 10) A,a,/J (n - r,| , z', z) (0| df 

Where Q a p is a quadrupole transition operator; the elec- 
tric field operator Fp is due to the Frenkel exciton as in 
(|11[> : Di_ at p is the Green's function in momentum space. 
The interaction parameter is given by the following tran- 
sition matrix element: 



T k = (W,-k\H int \F,k) = 
dzdru (W,k\H mt \F,k) 



EE 

a,[},i n 



The quadrupole transition matrix element may be 
obtained by using the relation ^ J2 <pis (q) = 

q 

fis ( r e — fft = 0) between the wave function of the rela- 
tive electron-hole motion of the IS* exciton ipis normal- 
ized to the unit are S and its momentum representation 
Vis.q- Hence, the Fourier component of the quadrupole 
exciton transition from the ground state of the crystal^ 
is given as: 



(W,\c\Q a , p \0) 



ikrii 



E^ 1S (q) 



2 XQ a p e 



/kr. 



7T Qb y/S 



Xl (z)x h k (z) 



Considering the weak dependence of the conduction u c , q 
and valence u v ,q Bloch functions on the wave vector q, 
the quadrupole moment of the inter-band transition may 
be written as: 

Q a ,0 ~ Q a ,p (q = 0) = — / dsn* (s) x a xpu v<i (s) 
vo J 



with the integration taken over the unit cell of the 
cuprous oxide vq. The expectation value of the FE exci- 
ton dipole transition operator is given as: 



(0\df (n)|F,k) = 



-a F dfe~^ 
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The Green's function in momentum space for the given 
geometry is given by 

A.c.,/3 (k, z', z) = -^L-k 2 e k ( z - z ') x 

„(* + «,)(S. + ^)(* i + v .) 

For the selection rules based upon the specific geometry 
of the problem and the corresponding Green's function 
see Table III. Without loss of generality we may take 
the wave vector along the x - axis. From the Green's 
function; polarization selection rules^; our specific type 
of interaction energy and the requirement of having real 
eigenvalues one can conclude that the only possible di- 
rection of the wave vector and dipole orientation is d% 
and Q x .z respectively. 

We note here that there is another possible approach 
to the problem of resonant interaction between these two 
types of exciton. If one modifies results of MoskalcnkoA2. 
for our case and introduce an effective polarization due 
to quadrupole transitions via a k dependent inter-band 
dipole element 

d y -,k ~ \Qij\ {eikj + ejki) (16) 

Then one uses the corresponding interband polarization, 
this would coincide with a case of 'dipolc-dipole' transi- 
tions, when the electric field created by the set of WE 
dipoles will interact with the dipoles in the organic. In 
this way one could directly use the results of Agranovich 2 ". 
Due to an additional requirement for the Hamiltonian to 
be hermitian, it would give us the following result for the 
coupling constant: 

(k) --. Re (d£ k D Ja d£) ~ \Q x ,j\4k xej e x 

Here we used the fact that the wave vector has only an x 
component. And so this coupling parameter vanishes in 
this approximation as j ^ x. For that reason in this arti- 
cle we study interaction of the cuprous oxide quadrupole 
with the gradient of the electric field created by organic 
Frenkel exciton placed near the interface between the or- 
ganic PS : CA and inorganic C112O. 

Calculating the integral of the coupling parameter and 
summing over the dipoles for the case when a perturba- 
tion approach is applicable yields: 

r (k) = Ti (k) + r 2 (k) 

Here Ti (k)- corresponds to the unperturbed coupling pa- 
rameter, and the last term is due to the perturbation. 

r x (k) - 8^ ke~ kz 'sinh(^) Q xz d F z 
1 6 (e (k) + i) L w ^ 1+ (fc£™) 2 ^ a F a b L w 

(17) 

To increase this coupling term one may consider the 
organic host PS:CA with DCM2 organized in a multi- 
layered structure. Generalizing the same approach to 



this case of multi-layered organic with distance between 
the layers to be equal to twice the radius of the Frenkel 
exciton (compact composition) it can be shown that one 

has to multiply (JT7J) by | 1 . The factor of two 

in the last expression indicates that the organic is placed 
on both sides of our IQW and one must consider the 
symmetrical Green's function in the interaction. 

Another term in the coupling parameter is due to in- 
duced Stark effect and has the form: 

512 LjkF k -C-M e .cosh(^) 
21 >~ 27 ^h 2 k^Li + 10n 2 k 2 L 2 w + 9^ 

Where Fk - is electric field induced by the organic layer 
of dipoles and the terms proportional to lP w have been 
omitted. 

For the region of wave vector and width of the IQW 
where one can not use perturbation theory another form 
of the coupling parameter can be introduced. In this 
case it is convenient to use the effective parabolic poten- 
tial and the envelope functions in the form (|15| . and the 
coupling parameter has the form: 

r(fc)~fc 2 e^'^exp(-i%p)x 

x (erf (2tt + kR )) + erf ( 2?r - kR ))) 

(18) 

Here for simplicity we put z e w — Zh, or in another words 
m e ~ rrih- For numerical estimation of the quadrupole 
matrix elements we used the well known^ result for cor- 
responding oscillator strength: 

3e^n \ZAa/ 

Where A represents affect of quantum confinement, and 
|k | = 2.62- 10 5 cm- x . 

The dispersion of the hybrid state in case of negligi- 
bly small detuning of IS exciton with the DCM2 tran- 
sitions (resonance) is presented on FigJS] The first (1) 
dashed lines correspond to the coupling parameter when 
the Stark effect is neglected. So it corresponds to to 
the upper part of the dispersion in this approximation. 
To show the effect of the 'parabolic' approximation on 
the coupling we drew the second dashed line (2). And 
the third dashed line is when the Stark effect is small 
enough to be taken as a perturbation without applying 
the 'parabolic' approximation. 

Note that we consider the quadrupole-dipole interac- 
tion only between Cu-iO and DCM2 excitons and ne- 
glect the coupling with PS:CA host transitions, as it 
has much smaller oscillator strength of the correspond- 
ing inter-band dipole moments for the transitions are off 
resonance with the inorganic. Neglecting the Stark effect 
from this numerical estimation it is easy to see that for 
the values of kR « 1 the expressions (fT8| and (fT7|) are 
equivalent. 
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FIG. 5: The solid lines represent upper and lower branches 
of the quadrupole-dipole hybrid dispersion when the coupling 
is calculated in the parabolic approximation and the induced 
Stark effect is taken into account; the dash lines correspond 
to the coupling parameter for different approximations: (1) - 
infinite IQW and the induced Stark effect is neglected, (2) - 
parabolic approximation, (3) - infinite IQW with the induced 
Stark effect treated as perturbation 

VII. RESULTS AND DISCUSSION 

Although the oscillator strength of the IS quadrupolc 
exciton is three orders of magnitude smaller than 
for the nearest dipole allowed exciton, strong spatial 
dispersion^ makes the maximum of the coupling param- 
eter comparable with those for the dipole-dipole case. 

The IS" quadrupole exciton has a rather big mass of the 
(w 3mo) and so is comparable to the mass of the Frcnkcl 
exciton (> 5too)- Also one has to take into account a 
significant effect of IQW width fluctuation. Indeed, the 
effect of the width fluctuation is energy change between 
confined and the bulk cases, which corresponds to energy 

drop Of Ei S ,binding (A = 1) - EiS,binding (A = 0.881) = 

17 meV (See FigJ3]). Hence , due to both of this effects we 
neglected the kinetic energy of the confined quadrupole 
15 exciton (the motion occurs in Frcnkcl like way, by hop- 
ping between sites of localization). Thus the quadrupole- 
dipole hybridization effect is more pronounced than the 
dipole-dipole hybridization because it is not masked by 
kinetic energy of the exciton. 

Also the hybridization effect is not masked by the large 
radiative decay rate of the Frenkel exciton. Instead of the 
classical spontaneous radiative rate of the organic we use 
the 'solid state solvation' process to obtain a different 
relaxation mechanism. Namely, one has the dynamical 
red shift of the DCM2 Frankel exciton during the 'slow' 
phase of the solvation process, when the energy of the 
photon excited DCM2 molecules is partially transfered 
to the polar CA molecules by non-radiative dipole-dipole 



cm. interaction followed by immediate (compared to the 
ts) radiative decay. To compare the radiative decay rate 
(life time) with the hybridization parameter we simplify 
the dynamics to the following statement. The Frenkel ex- 
citon has fixed resonant energy with the IS quadrupole 
Wannier exciton but has effective lifetime equal to rs>2£. 
So if one omits all the possible non-radiative channels of 
decoherence except the radiative decay than the radia- 
tive decay rate from the hybrid to the ground state is 
calculated using Fermi's golden rule: 

frlhybrid (u, I; k) = 2ir \ (g\ H mt \u, I; k)\ 2 5 (hcu k - E Ut i) = 

A fu h B lu n hT s +T 1S 
T s T\S 2 T s Tl S 

From this estimation the ratio of the hybrid radiative 
rate to the maximum of the coupling parameter is 0.0006 
which is two orders of magnitude smaller then predicted^ 
value 0.09 for the dipole-dipole hybrid in a quantum 
do1j2£. Indeed it was shown in case of fully localized 
quantum dot dipole-dipole hybridization that the ratio 
of the hybrid exciton damping parameter to the coupling 
is proportional to lh ""^ nd oc {y/ab - in case of quan- 
tum wires^) so quadrupole-dipole hybridization takes ad- 
vantage of the small radius of the IS quadrupole. Also 
the fact of comparable life time of the organic and 15* 
quadrupole exciton significantly influences the ratio. 

A noteworthy feature of the resulting upper (u) and 
lower (1) branch dispersion in Fig. [5] is that the well pro- 
nounced minimum on the lower branch can be populated 
by pump-probe experiment and now it may be possible 
to have finite critical temperature in case of quasi 2D 
excitons, due to the fact of non parabolic dispersion of 
the lower branch of the hybrid. This can provide a good 
basis for searching for BEC in such a hybrid. A new type 
of an interface polariton may also be expected. 

Note that we mentioned the possible BEC as a possi- 
bility only. Judging from the fact of the big saturation 
density of the hybrid and minimum of the energy on the 
lower branch, this minimum would make the hybrid ex- 
citon strongly localized and many dissipative processes 
(such as Auger heating^) are going to be suppressed. 
Also as we mentioned the non-parabolic dispersion al- 
lows the system to have a finite condensation tempera- 
ture contrary to the case of 2D excitons with parabolic 
dispersion. But the main question about the life time of 
the hybrid is still a subject of our current research. So 
we do not make a concrete statement about the BEC of 
the hybrid but rather make an educated guess on that. 
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TABLE I: Selection rules for the excitons in bulk Cu 2 0. F-forbidden; A-allowed transition 
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TABLE II: Selection rules for the exciton in a C112O quantum well. F-forbidden; A-allowed transition 





Qxy 


Qxz 


Qyz 


k [i,o,o] 


Dxxy — 
■Dyxy — 
Dzxy 


D xxz = _^fc a e fc l«-'),d = (d„0,0) 

DyxZ 

D zxz =i£ I k 2 e k ( z - z '),d= (0,0, d.) 


D X yz 

Dyyz 

D zyz = 


k [o,i,o] 


D X xy — 
■Dyxy — 
Dzxy 


Dxxz — 

DyxZ 
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D X yz — 

Dyyz = -^k 2 e k ( Z - Z '),d = (0,^,0) 

Dzyz = i^k 2 e k ( z - z ') , d = (0, 0, d z ) 
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D X yz — 
Dyyz = 
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TABLE III: Green's function selection rules 



